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Abstract 

Using a generalized transfer matrix method we exactly solve the Schrodinger equation in a 
time periodic potential, with discretized Euclidean space-time. The ground state wave function 
propagates in space and time with an oscillating soliton-like wave packet and the wave front is 
wedge shaped. In a statistical mechanics framework our solution represents the partition sum 
of a directed polymer subjected to a potential layer with alternating (attractive and repulsive) 
pinning centers. 
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Schrodinger equation plays a central role in modern physics. In the Euclidean space-time, 
the time is imaginary and the equation has interpretation in statistical mechanics. In contrast to 
the stationary case, i.e. equations in time independent potentials, where many exact solutions 
exist, for the general time- dependent situation exact results have been lacking so far. 

In this work we present an exact solution, albeit in the Euclidean lattice space-time. We 
consider a periodic potential alternating between attractive and repulsive, placed at the space 
origin x = 0. The quantum analogue would be that of a particle in a forced potential with 
energy pumped in and out periodically. This problem has several links with important branches 
of modern theoretical physics. In the statistical mechanics framework our system, with a delta- 
like potential, is usually employed in the lattice models for wetting and depinning transitions of 
directed polymers jjj], as well as it represents a simplified version (but nevertheless non trivial) of 
the KPZ equation for kinetic interface growth 0. In a pure QM context the present problem has 
applications in: wave function collapse in mesoscopic systems, Mott-type hopping, transport in 
disordered systems, quantum electronics, time-dependent Zener tunneling (see || and references 
therein). Moreover it has been recently pointed out that the Schrodinger and diffusion equations 
are very deeply linked if represented on a lattice. In fact one can prove that, apart from the 
usual analytic continuation, they represent two aspects of the same probabilistic problem on 
a lattice manifold Q. This maybe provides a new interesting interpretation of non-relativistic 
QM on a lattice. 

The strategy of our approach is the following: we assume that at large times the wave 
function has the same periodicity of the potential. We also require that it is normalized for all 
times. For the symmetric case, i.e. when attractive and repulsive potentials have equal strength, 
the wave function is localized around the origin ||. The normalized solution oscillates in the 
time direction and its amplitude is weaker and weaker when the spatial distance from the origin 
is larger and larger. Our wave function also shows that the wave front has a wedge shape as 
would be required by causality: the potential influence can only travel with finite velocity. 
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Let us start by considering the imaginary-time <i-dimensional Schrodinger equation 

d t ^(x, t) = (V 2 + 5{x) cos(wt)) ¥(a?, t), (1) 

with the normalization J \*f>(x, t)\ 2 dx = 1. Here we have considered the symmetric case: the 
time average of the potential is zero. However our general solution also includes the nonsym- 
metrical situation. For sake of simplicity we will use the 1-d notation for variables and operators 
even if we can work in a high dimensional space. 

It can be shown that the above equation can be derived from a Hamiltonian ^,|| which 
describes an elastic chain, or a directed polymer, in a periodic potential (alternatively attractive 
and repulsive). The Hamiltonian reads 

L L 

n L {{h k }) = JJ2 \hk+i ~ h k \ p - u k$h h ,o, (2) 
k=l k=l 

and, more specifically, it gives the energy of a polymer in a (d + l)-dimensional space under 
the competing effect of the pure tension term proportional to J and the potential layer at the 
origin. Here h k defines the position of the polymer at "time" k, while u\~ is set to u > for 
even k and — v < for odd k. In the usual lattice version, one introduces RSOS conditions, 
that is the height difference — h k \ can only take values or 1 and overhangs are forbidden. 
Mathematically speaking, this is completely equivalent to assume p = oo [Q] . The first term in 
(1) discourages large humps of the polymer (it corresponds to the laplacian in (1)), while the 
second term gives a positive (resp. negative) contribution to the total energy when the particle 
passes through the origin for odd (resp. even) times k. Hence it is the discretized (and non 
symmetric) form of the continuous sinusoidal potential in the original Schrodinger equation. 

The equilibrium state of the polymer, whose energy is defined in (2) for a given configuration, 
is the result of the competition between two contrasting effects. At temperatures T > we know 
from general arguments [|J that large deviations from linearity increase the configurational 
entropy, but these jumps have also the effect of increasing the internal energy and a phase 
transition occurs at a given T c when the cost in energy is exactly balanced by the gain in 
entropy. This mechanism is completely specified by the free energy (per unit length) / = F/L. 
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From a geometrical point of view, below T c , the polymer is localized in the sense that its mean 
variance from the origin remains finite in the thermodynamic limit L —* oo. The opposite 
stands at T > T c in which the directed line behaves as a pure random walk in a d-dimensional 
space and performs thermal wandering. The QM counterpart of this phase transition is the 
binding/unbinding transition between eigenstates of the Hamiltonian defined in (1). After 
having clarified the mapping, our goal consists in the search of the partition function: 

Z L {x) = J2 e- Wi « fc *»/ T , (3) 
{M 

which is the discretized Feynman path integral representation of the wave function associated 
to (1). 

We then start by introducing the evolution equation for Zl(x) by means of a transfer matrix 
approach: 

Z L+1 (x) = [1 + (ol+i - 1)<5 Xi0 ] (4) 

L 

x Z L (x) + t^(Z L (x + l)+Z L (x-l)) . 

i=l 

The search of the ground state of eq.(l) is equivalent to the determination of the stationary 
solution of this recursion equation for L — > oo. In the thermodynamic limit we know that only 
the largest eigenvalue e of the transfer matrix (4) gives a significant contribution to the free 
energy density, that is we have f — — loge . We recall that, due to the particular form of the 
potential layer, we expect to find two distinct ground state wave functions in the bound regime, 
one for even and one for odd times. 

The first step consists in the search of a self-consistent expression for the partition "wave" 
function. The full calculation has been already carried out in a general hypercubic lattice ||: 
the result is that for any set of free parameters {d, u,v,J} the normalized partition functions 
in the Fourier space Gf(k) = Zf(k)/{l + 2dt) L (with t = exp(-J/T)) @ satisfies the following 
implicit equation 

Gt +1 (k) = C%t(k)+A f Q d d qC%(q)t(q), (5) 



where A = exp(u/T) — 1 = A + for even times and A = exp(— v/T) — 1 = A~ for odd times. 
Moreover we have introduced the quantity = [1 + 2t J2i=i cos(7rA;i)]/(l + 2dt). To take 
into account the periodicity of the potential layer we apply transformation (5) twice, so that 
we can relate G , ^ i+2 to G 2L and G 2L+l to G^^- By using the ansatz, valid at large L, 
G±(k) = e 2 G±_ 2 (k) we finally get § 



+ -f^K*. (6) 



where we have introduced the two constants 

K+= f 1 d d qC(q)G + (q), n = 1,2. (7) 
Jo 

We simply find, by symmetry, that G-{k) can be derived from G+(/c) by interchange between 

A + and A~ . A phase transition occurs when the maximum eigenvalue converges towards 1 + 

H , and it can be studied by means of the mass gap fi defined as the inverse of the transversal 

correlation length, i.e. ^(x) ~ exp(— fj,x). One can prove that e ~ 1 + fi 2 near the transition 

point H . The r.h.s. of eq.(6) implicitly depends on G ± (k), the full knowledge of the wave 

functions would therefore require 4 independent equations for the constants K^ 2 - The first 

couple of equations can be simply derived from (6) by multiplying both sides for £(/c) and 

integrating over the fc-momentum. At the end we obtain the following identity 

/ A ± A^ \ 

Kf A^l 2 + — - -2i - 1 + KfA+l! = 0. (8) 



where 



l + 2dt 



Ii = ^(/i^-flifO)), 



12 = ^^(/i(0) + 3i(0))-l, (9) 



and 



poo 

f e (x) = / due-^ u I {xl (u) d , 
Jo 

poo 

g £ {x) = / due- £ * u lU-u) d . (10) 
Jo 
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In the above we have introduced the two quantities £\ = [s(l+2dt) — i\/2t, £2 = [e(l+2di) + l]/2t 

and we have used an integral representation of the modified Bessel function of integer order 

In{v) = f q ex.p{u cos y} cos ny dy . It is simple to convince oneself that any other attempt to 

find close equations for these constants from (6) would lead to relations which are not linearly 

independent respect to (8). 

For the meanwhile we neglect the search of the second couple of equations and we concentrate 

on equation (6), which must be antitransformed in order to get the real space form of our 

solution. In performing the calculation we encounter non-trivial integrals of the form: 

ri d f(k) n 
J n (x)= d d kT\cos(TrkiXi) 2 V v ' » n = 0,l. (11) 

Without entering into mathematical details, we simply note that, by means of an appropriate 
Feynman integral representation, we can simplify them and after some more algebra we finally 
obtain that 

Mx) = [/«(*) +*(*)], 

Jl{x) = i±M [/£ ( x )_ 5e(x)] . (12) 

The symmetry of the system tells us that any function of x must be invariant under axis 
reflection; this is why in the above integrals only the absolute value of x appears. Note that 
after this manipulation, the dependence of our wave function on the spatial variable x comes 
only from the order of the Bessel function involved in the integrals. Moreover, by using the 
property that I\ x \(u) is an even (odd) function of the argument for even (odd) values of \x\ it 
is clear that our ground state solution has a particular oscillation as a function of \x\. As we 
will discuss below, this is nothing but a direct consequence of the alternating potential we have 
put at the origin. We also note that f £ (x) is a well known integral in the statistical mechanics 
context H : if one takes e = 1 (or £\ = d) , it is associated to the total probability that a random 
walker, started from the origin, could finally reach a point a; in a d-dimensional cubic lattice. 
The divergence of fi(x) for d < 3 has also important consequences in the polymer depinning 
framework ||. 
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To summarize, our non-normalized wave function reads 

¥±(*) = {^Kf + ^^Kpj Jo(x) + A+KfMx). (13) 

The above formula gives the qualitative behavior of our solution, but to get the full normalized 
wave function we also need a closed form for K^ 2 - As we know that e > 1 for a non zero 
mass gap, we can restrict to the localized phase and impose the correct normalization of the 
wave function. Therefore we ask that J2{x}^±( x ) 2 = 1 in the real space, or equivalently, 
Jq d d kG±(k) 2 = 2 in the momentum space. This condition would give us the second couple 
of equations we needed to find Kf 2 . Unfortunately the integrals resulting from (6) after this 
manipulation are much harder to handle with. In fact we find that the normalization condition 
reads 

F 2 (Kf ,A^)M 2 + F 1 (Kf 2 ,A ± ,A^,d, t)M t + 

F (K^ 2 ,A ± ,A^,d,t)M = 2, (14) 

where Fq,Fi,F 2 are some algebraic functions of their arguments and 

M n = C d d k 9 , n = 0,l,2. (15) 

We believe that no simple way can be found to simplify these high-dimensional integrals, and 
then at this point we are forced to restrict ourselves to the l-d case. In a pure statistical mechan- 
ics framework one could, in principle, ask ^±(x) to be sommable and not square-sommable. 
In this case we can perform all integrations and get the complete solution in all dimensions. 
Nevertheless, in the spirit of QM, in the present study we prefer to use the usual normalization 
of the square of ^±(x) . At d=l all previous integrals can be exactly solved and we find that 

fe(x) = 



9e(x) = j ) —. (1(5) 

e 2 2 -l(e 2 + Je 2 2 -l 



From these formulas we get immediately the result that the mass gap in the localized phase is 



given by ll = log(ei + \J e\ — 1) and then it vanishes, near the transition as y/e — 1. This is, a 
posteriori, the proof that the limit e — > 1 + gives the transition point, as previously stated. 
By using (8) and (14) we can finally solve for Kf 2 '- 



Kf = V21 l {[l-2A^l 2 + (A^fll 



M 



+ 



2A T 1 1 - 2(A T ) 2 1 1 1 2 ] Mi + (A T ) 2 lfM 2 ] 



-1/2 



* 2 ~ {A^ A±T X l + 2dt) Kl - {7) 

We now have the complete solution of our problem, since the ground state (bound) wave function 
is given by (13) with the coefficients defined in (17). In particular we find that, at d=l, 

M 2 = ■^(h 1 + h 2 + e(h+l 2 )), M 1 = —(h-l 2 ), 

M = + Z 2 )-/ ll -h 2 ), (18) 

with 

hi* = ^(4 2 -ir 1/2 , 

h,2 = (l^ eii2 (4 2 -l)-3/=>. (19) 

After rearrangement of all the quantities we can now have a look of our result. The 1-d wave 
functions ^+(x) and \V-(x) are plotted in Fig.l for a given set of parameters {u,v,t} in the 
localized phase. In order to check out the validity of our solution, we have performed some 
numerical simulations, and for any set of parameters we have found perfect agreement with the 
above analytical solution. 

The information of the alternating potential at the origin propagates in the transversal 
direction with finite velocity, in a solitondike fashion. This means that if we look at the shape 
of our solution in the space-time manifold we expect to find that the oscillations due to the 
alternating perturbation of the potential decrease by increasing the distance \x\ from the origin 
and the wave front has the shape of a wedge (see Fig.2). 



Performing the limit e — ► 1 + we attain an unbound state: the peaks at the origin disappear 
and the two wave functions merge one over the other at diverging transversal correlation length 
(i.e. at vanishing = If one is interested at the behavior of the "gap" Ao = l I'+(0) — 

^_(0) near the transition to an unbound state, then all information comes from the exponent 
a, defined as Ao — (e — l) a for e — > 1 + . From (13) and (17) we simply find that for (s-1) < 1 
the gap reads Ao ^ + 2t) / A(A + — A~K^~) and then two situations are possible: (i) 
we can converge to the critical point for generic A + ^ A~ or (ii) we can attain a delocalized 
phase for A + = A" . More precisely, in the second hypothesis, we approach the unbound state 
on the manifold u = v and u — > linearly with e — 1. This difference is not trivial but rather 
reflects a physical property of our system: at d = 1 it has been proved that in the two above 
cases the phase transition is of 2nd. and 4th. order, respectively ||. By taking the dominant 
contributions of the integrals defined in (9), (12) and (18) for e —* 1 + in the formulas (17) for 
K j , after some calculations we find that Kf ~ (e - l) 1 / 4 . This result is the same in both cases 
described above, since the limit u = v — > does not modify the leading order contribution in 
(17). Then we finally obtain that, depending on the way we approach the critical state, a = 1/4 
(i) or a = 5/4 (ii) (recall that A ^ linearly converge to for e — > 1 + in the last case). We finally 
remark that our solution, in a pure statistical mechanics context, has a certain relevance because 
it is the partition function associated to Hamiltonian (2), which describes, as above explained, 
the energy of a directed polymer in a periodic potential, a problem with several applications in 
statistical physics 

In conclusion we have solved an imaginary time Schrodinger equation with a non-trivial 
time dependent potential on a lattice. The critical behavior separating localized and delocalized 
phases has interesting properties characterized by non trivial scaling exponents. We expect that 
the counterpart in continuum would have qualitatively same behavior, which may have in turn 
wider implications in the non-equilibrium quantum physics. 
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Figure captions 

Fig. 1 

The wave functions ^>(x) = and $(x) = versus x for the 1-d case in the bound 

state. The curves are obtained by smooth interpolation among the points given by the analytical 
solution on the lattice. Numerical simulations fit exactly the above curves for any set {u, v, J} 
of free parameters. 

Fig. 2 

The full space-time shape of our wave function t). As explained in the text, the wave front 
is wedge-shaped. The two functions (resp. are easily found by intersecting the surface 
with a t = const, plane for times t such that ^(x = 0, t) reaches a peak (resp. a valley). 
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